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We derive a quantum master equation which describes the dynamics of the ensemble-averaged 
state of homogeneous disorder models at short times, and mediates a transition from coherent 
superpositions into classical mixtures. While each single realization follows unitary dynamics, this 
decoherence-like behavior arises as a consequence of the ensemble average. The master equation 
manifestly reflects the translational invariance of the disorder correlations and allows us to relate 
the disorder-induced dynamics to a collisional decoherence process, where the disorder correlations 
determine the spatial decay of coherences. We apply our theory to the (one-dimensional) Anderson 
model. 
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I. INTRODUCTION 

It was the insight of Anderson that disorder can sub¬ 
stantially modify the dynamical behavior of quantum 
particles: The destructive quantum interference due to 
multiple scattering off impurities in the wire potentially 
brings the electrons to a halt, giving then rise to An¬ 
derson localization [T]. Even when the consequences of 
disorder are less drastic, its interplay with quantum inter¬ 
ference can still alter the mobility pattern, causing, e.g., a 
transition from ballistic propagation to weak localization 
[2]. While these interference effects already occur on the 
level of single disorder realizations, they even prevail un¬ 
der an average over many disorder realizations, this way 
stripping off individual peculiarities and defining a statis¬ 
tically robust effect. Anderson localization, e.g., unveils 
its characteristic trait, exponential wave function tails, 
on the level of the ensemble average. 

The possibility to implement disorder models with 
highly-controllable cold atomic gases has brought it into 
reach to access disorder phenomena and their underly¬ 
ing quantum origin even on the level of the spatially re¬ 
solved atomic density n{r,t) [31 H]. It was, for instance, 
observed that the ensemble-averaged correlation function 
of density fluctuations exhibits, at long times, character¬ 
istic long-range correlations, which can be traced back to 
the macroscopic coherence in the gas 1313]. Here, we in¬ 
vestigate the evolution of quantum coherence under the 
disorder average at short times. We find that the spatial 
pattern of the coherence loss of the ensemble-averaged 
state is directly related to the correlations in the disor¬ 
der potential. This loss already happens at ballistic times 
much shorter than the mean free time r, where the dis¬ 
order does not yet have a significant effect on the level of 
single realizations. 

We emphasize that this effective decoherence of the 
ensemble average state does not correspond to a loss of 
information as it generally occurs in the presence of an 
environment. In our case, single disorder realizations fol¬ 
low the unitary dynamics of isolated quantum systems, 
i.e. the occurrence of quantum interference phenomena 


which survive the ensemble average, such as Anderson 
or weak localization, remains untouched. The coherent 
nature of the dynamics of single realizations can for in¬ 
stance be recovered by considering higher-order corre¬ 
lators, such as, e.g., intensity correlations. The loss of 
coherence of the ensemble-averaged state, on the other 
hand, is a consequence of the fact that different disor¬ 
der realizations propagate an initially pure state into dif¬ 
ferent evolved states, and that their averaging generally 
results in a mixed state. 

To establish our results, we derive a general Lind- 
blad master equation for the evolution of the disorder- 
averaged state on short time scales, allowing us to inves¬ 
tigate the transient dynamics for arbitrary initial states 
Pq. In this approach, the dynamical impact of the disor¬ 
der is reflected by the structure of the resulting master 
equation. In particular, coherent and incoherent contri¬ 
butions to the dynamics of the ensemble-averaged state 
are consistently separated. As we show, the evolution 
generated by the master equation for the one-dimensional 
(ID) Anderson model perfectly agrees with the short- 
time dynamics of numerically exact simulations thereof. 

Let us stress that our approach lies at the interface 
between quantum transport theory of disordered sys¬ 
tems and the theory of open systems. It complements 
other perturbative methods to treat disorder dynamics, 
e.g. based on averaged propagators and/or diagrammatic 
methods IZIIE]. Alternative evolution equations for the 
ensemble average state have been proposed in Ellin]. 

A comprehensive understanding of the disorder- 
induced dynamics at short times is also of practical rel¬ 
evance, as it permits one to access the detrimental im¬ 
pact of perturbations on the functioning of quantum de¬ 
vices. To see this, let us consider a simple example, the 
double-slit experiment. There, the observed fringe pat¬ 
tern, which represents the purpose of the device, strongly 
relies on the delicate interplay between the delocalized 
state prepared by the slits and the phases accumulated on 
the way to the screen. What happens if the particles are 
disturbed along their way, e.g. if they propagate across 
a disordered scattering potential towards the screen? As 
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we show in Fig. 1, averaging over many realizations of 
the disorder potential gives rise to a continuous-in-time 
decay of coherences, i.e. the visibility of the interference 
pattern in momentum is monotonously reduced as time 
elapses. In other words, while single realizations exhibit 
distorted interference fringes, the ensemble average re¬ 
covers the structure of the undisturbed pattern, however 
with an increasingly reduced visibility. 
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FIG. 1. (Color online) Decoherence dynamics induced by the 
disorder average: Evolution of an initial spatial superposition 
state in the (one-dimensional) Anderson model, mimicking 
the double-slit experiment in the presence of disorder, (a) 
Density matrix of the spatial superposition of two Gaussian 
wave packets a.tt — Oh/J (initial state), (b) Ensemble average 
state Pens(t) a,t t — 0.8h/J (J denotes the hopping constant), 
with disorder strength W = 5J, and averaged over K = 100 
realizations. One observes a decay of the coherences between 
the two peaks, as well as of each individual peak. The loss of 
coherence is also reflected in the reduced visibility of the in¬ 
terference pattern in momentum (c - blue solid a,tt — 0.8h/ J, 
orange dashed at t = Oh/J), and in the decay of the purity 
Pens(t) = tr[pens(t)^] of the State (d). 


II. HOMOGENEOUS DISORDER 

We consider a single quantum particle subject to a 
homogeneous disorder potential, i.e. correlations among 
different locations are translationally invariant. For sim¬ 
plicity, we focus here on a one-dimensional, discrete (in¬ 
finitely extended) configuration space, comprised of sites 
|j) with lattice spacing a; however, as will become clear 
in the course of the article, our theory works as well 
for continuous configuration spaces, higher dimensions, 
finite-size topologies, and many particles. In the case 
considered here, the Hamiltonian for a single disorder re¬ 


alization is given by 

+1)01)+E"0j)OI: (1) 

jez jez 

where the (infinite-dimensional) vector e comprised of the 
random on-site energies £j distinguishes different disor¬ 
der realizations. The tunneling or hopping term is char¬ 
acterized by the tunneling/hopping constant J, which 
controls the maximal propagation speed in the system. 

The distribution p(e) of the on-site energies is as¬ 
sumed to be homogeneous. Besides the normalization 
{riigz/ = I, we thus require that the expec¬ 

tation values and two-point correlation functions satisfy 



Note that we assume without loss of generality that the 
expectation value e of the onsite energies vanishes. For 
convenience (as will become clear below), the two-point 
correlation function C'(Aj) is measured in units of the 
hopping constant J (as well as all other quantities with 
the dimension of an energy). As we will see, it is sufficient 
to characterize the expectation values and the two-point 
correlations of the disorder distribution, as only these ap¬ 
pear in the master-equation description at short times. 
Of course, homogeneity requires all higher-order correla¬ 
tion functions to be translationally invariant as well. We 
restrict ourselves to disorder which is diagonal in the site 
basis; however, other forms of homogeneous disorder are 
also conceivable. 

In the case of the Anderson model [T], the on-site 
energies of different sites are completely uncorrelated, 
i.e. the disorder distribution p(e) decomposes into a 
product of identical single-site distributions, p(e) = 
Yii^zPsi^i)- The box-shaped single-site distributions 
Psis) = 0(IF/2 + e)0(IF/2 — e)IW are characterized by 
the disorder strength W (0 denotes the Heaviside func¬ 
tion). The translation-invariant correlation function is 
accordingly given by 
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(3) 


III. SHORT-TIME EVOLUTION. 

We now derive a quantum master equation which ac¬ 
curately describes the ensemble average dynamics of dis¬ 
order models such as Q at short times. In particular, it 
renders the, in general, incoherent nature of the ensemble 
average dynamics manifest in terms of the emerging Lind- 
blad terms. To this end, we first consider the unitary time 
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evolution for a single realization of the disorder potential, 
pg{t) = Ug'{t)poUl{t), with the initial state po at tg = 0 
and the time-evolution operator U^{t) = exp{—iH^t/h). 
Since we are interested in the evolution on short time 
scales, we expand to second order in dt: 

Pe-(dt) =Po + [po, He] (4) 

+ ^ [sePoHe- Injpo - IpoHI^ + Oidt^). 

The second-order term is structurewise reminiscent of a 
Lindblad term, and, indeed, upon averaging over differ¬ 
ent realizations, the leading incoherent contributions to 
the disorder dynamics arise at second order in time. To 
see this, we take the average p — {rijgz I dci} of 

which yields 


p{dt) = Po -k ^dt [po, .ff] -k dt I J de, | 


(5) 


1 


1 


X HePoHe- -Hipo - -poHi + 0{dt^). 


In the first-order von Neumann term we exploited that 
the initial state is independent of the disorder re¬ 
alization. It therefore commutes with the ensemble 
average, resulting in the average Hamiltonian H = 
{rizGZ J dsi} p{^Hg-. Such reduction is, in general, im¬ 
possible for the second order term, which ultimately gives 
rise to incoherent dynamics. 

To convert ([^ into a differential equation of Lindblad 
form, we must not restrict our treatment to the leading 
contributions in dt, since we would thus lose the inco¬ 
herent part of the dynamics and end with the coherent 

evolution induced by the average Hamiltonian H alone. 
To consistently identify next-to-leading order contribu¬ 
tions, we replace Hg- —)■ {Hg— H) + H. Equation ^ can 
then be rewritten as 

p(dt) =po -f ^dt [po, H] (6) 

+ —i^HpoH--H po-^PoH j 

+ I n / I - ^)PoiHe - Ti) 

- +0(dt3), 

where the first two lines represent the von Neumann com¬ 
mutator, and the last two lines the Lindblad terms, re¬ 
spectively, each to second order in time. The decoherence 
rates associated with the Lindblad terms increase linearly 
in time. 

It follows that Eq. ^ solves, to second order in time, a 


Lindblad master equation for the ensemble average state, 

^ = - ^[H,p\ + I J de,| 7r(t) (7) 

X (^LgpLi- ^LiLgp-^pLiLe^ , 

which captures the disorder dynamics at short times. 
The (time-independent) Lindblad operators Lg and the 
corresponding (time-dependent) decoherence rates ‘jg{t) 
read 


Lg = 


Hg-H 

Eo 


le{t) = 




( 8 ) 


where the characteristic energy scale Eq is introduced 
in order to obtain the appropriate dimensions; as stated 
before, in the case of the model Q it is conveniently cho¬ 
sen to be the hopping constant J. We thus find that the 
ensemble average accounts for each disorder realization 
by an independent Lindblad term, where the Hermitian 
Lindblad operators are given by the offset of the disorder 
Hamiltonian from the average Hamiltonian. The asso¬ 
ciated decoherence rates are proportional to the proba¬ 
bility p{e) for the realization to occur and scale linearly 
in time, i.e. the rates vanish at t = 0. The latter ex¬ 
presses that there is no incoherent contribution to the 
dynamics at first order in time. The validity range of 
the short-time approximation Q depends on the com¬ 
position of the underlying disorder ensemble and must 
be determined case by case. While the master equation 
0 does not require, e.g., weak disorder, the time scale 
on which it yields reliable predictions in general depends 
on the disorder strength. Below we will give a numerical 
estimate for the Anderson model. 

We emphasize that the disorder master equation Q 
still holds for arbitrary systems and general disorder dis¬ 
tributions, since we have not yet made use of the Hamil¬ 
tonian 0 and/or of the homogeneous distribution 0 - 
In the Appendix, we thus evaluate the short-time disor¬ 
der dynamics 0 for two unrelated, yet instructive exam¬ 
ples: a particle of mass m in one-dimensional, continuous 
space, subject to a random i) linear or ii) harmonic poten¬ 
tial. In these cases one hnds that the short-time dynam¬ 
ics of the ensemble-averaged state are governed by the 
well-known Caldeira-Leggett master equation nn [la¬ 
in the case of the disorder model 0> the average 
Hamiltonian is given by the discrete hopping term, H = 
+ 1| + 1/ + l)(il) (we used Q), and the 
Lindblad operators are given by the disorder potentials, 
Li- = (with Eq = J). This is already 

conceptually appealing, since it demonstrates that the 
Lindblad operators are diagonal in the site basis; more¬ 
over, Eq. 0 predicts (confirmed by observation) an ini¬ 
tially quasi-free, dispersive evolution of the ensemble av¬ 
erage state in addition to the loss of coherence. However, 
this representation is not yet viable from a practical point 
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of view, in the sense that it is not amenable to transpar¬ 
ent approximations or efficient numerical simulation. In 
the following, we derive an alternative representation for 
homogeneous disorder models which resolves these is¬ 
sues and, in addition, reveals a connection to collisional 
decoherence. 


IV. COLLISIONAL DECOHERENCE MASTER 
EQUATION. 


behavior of (101 by neglecting the coherent dynamics ac¬ 
cording to the von Neumann commutator and solving the 
remaining equation in the position representation. One 

then obtains {j\pit)\f) = exp - j') \ (j|pol/), 


where the localization function F{j — f) (not to be con¬ 
fused with the exponential localization of the particle 
density in the Anderson model) follows from a Fourier 
(back-)transformation of G{q) and evaluates as 


To obtain an alternative representation for the short- 
time dynamics 0 of the homogeneous disorder models 
0 , we exploit that Lindblad master equations are invari¬ 
ant w.r.t unitary transformations of the Lindblad oper¬ 
ators. In our case this corresponds to a transformation 
from the position to the momentum basis. To this end, 
we perform the disorder integrals in 0 and are left with 
the double sum over the sites appearing in the Lindblad 
operators, 

-\\j){3\3'){j'\-p-\-p\3){3\3'){3'\), (9) 

where we already made use of the translational invari¬ 
ance 0. If we then rewrite the correlation function 
C{j) in terms of its Fourier transform G{q)^ G{j) = 
(q denotes the lattice spacing), we 
obtain with x = J 2 jez 30 -\j){j\ 


p = -iW.p] 


, 2JH 


hj2a 




dqG{q) (e*«"/^pe-‘«"/''-p) . 


— h}2a 

( 10 ) 

This is our main result. We find that the ensem¬ 
ble average dynamics of homogeneous disorder mod¬ 
els ([^ are at short times described by the discrete 
version of the collisional decoherence master equation. 
The (non-Hermitian, but unitary) Lindblad operators 
Lg = exp{iqx/h) describe momentum kicks, whose oc¬ 
currence is weighted by the momentum transfer distri¬ 
bution G{q). The latter follows by Fourier transform 
from the two-point disorder correlation function C(j), 
G{q) = {a/h)^^^^exp{-iqja/h)G{j) [H]. 

The collisional decoherence master equation (10) is 
usually known from an open-system context 
where it describes the decoherence that a heavy test par¬ 
ticle undergoes due to scattering in a background gas 
of light particles, i.e. no appreciable energy exchange 
occurs. It represents the simplest manifestation of a 
translational-covariant Lindblad master equation nsi- 

m. 

The master equation (10) allows us to deduce the de¬ 
coherence dynamics of the homogeneous disorder models 
([^. One can best understand the spatial decoherence 


/ /i/2a 

dqG{q)-GU-j'). (11) 

-hl2a 


We thus find that the disorder two-point correlations 
C{j — j') directly translate into the spatial decay of co¬ 
herences, in the sense that the stronger the correlation 
between two sites, the longer their coherence survives. 
This again reflects the fact that the disorder, i.e. the de¬ 
viations among different ensemble members, gives rise to 
the decoherence-like behavior. 


In the case of the ID Anderson model with the cor¬ 
relation function ([^, one obtains a constant momentum 
transfer distribution, G{q) = y^j 2 , and the localization 
function reads F{j — j') = ]^^(1 — be., while 

the populations remain unaffected, all spatial coherences 
undergo the same decay, independent of the separation of 
the two respective sites, since the sites are uncorrelated. 


As a second example we consider a Gaussian ran¬ 
dom potential with the correlation function G(j — j') = 

^ exp (— ° ), where ^ denotes the correlation 


strength and L the correlation length. Such correlations 
may, for example, emerge from a collection of Gaus¬ 
sian scattering potentials v{j — jn) with randomly dis¬ 
tributed scattering centers jn- This then yields the mo¬ 
mentum transfer distribution G{q) = exp 


and the localization function F(j — j') = ~ exp — 

^, i.e. there is a smooth, Gaussian transition 

into the regime of constant decoherence (|j — j'\a L). 
The coherence loss at short times caused by such Gaus¬ 
sian disorder correlations was also investigated in [20] 
in terms of path-integral techniques (for Gaussian initial 
states in the continuum and a harmonic average potential 



F[). In our language, the authors derive the localization 
function F(a;-a;') = (C/j2)(l-(l+2(3^)2)-i/2)^ which 
coincides in the short-range region (|a; — a;'| < L/2) with 
our result and shows qualitatively the same behavior in 
the long-range region. As the short-time master equation 
Q is derived without reference to a Hilbert space basis 
and therefore holds over the range of all sites, we inter¬ 
pret the quantitative deviation in the long-range region 
in terms of a breakdown of the path-integral approach. 
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V. NUMERICAL COMPARISON. 


In order to estimate the range of validity tmax of the 
short-time disorder master equation (10), we compare 


it in case of the ID Anderson model to the numerical 
ensemble average over a finite sample of disorder realiza¬ 
tions. In Fig. we show, in terms of an initial Gaussian 
state and for strong disorder with W = lOJ, that the 
master equation correctly predicts (relative error ±5%) 
the spatially homogeneous decay of the coherences up 
to about tniax = 0.2H/J, where the state has lost about 
45% of its initial purity p = tr[p^]. Similarly, one ob¬ 
tains for disorder strengths W = IJ and W = O.IJ va¬ 
lidity ranges of about tmax = 0.9Ii/J at a purity loss 
of 10% and tmax = Qfi/J at a purity loss of 1%, re¬ 
spectively. A more detailed analysis confirms that tmax 
roughly scales inversely with W, tmax oc 1/VF, or, in 
terms of the mean free path i, tmax oc Vi (similarly, the 
momentum-independent decoherence rate 'y{t) scales in¬ 
versely with the mean free time t, 'y{t) cx t/r). This 
suggests to interpret our theory in terms of an expan¬ 
sion in Wt. Notwithstanding, we can probe the regime 
of strong decoherence as induced by large W. Let us also 
emphasize that we could have chosen any initial state for 
this analysis. 


VI. EXPERIMENTAL VERIFICATION. 

Besides the conceptual insight provided by our theory 
into the incoherent ensemble average dynamics of dis¬ 
ordered quantum systems, direct experimental verihca- 
tions thereof are conceivable, for example based on ex¬ 
periments with ultracold atoms subject to optical speckle 
potentials. These systems have already been success¬ 
fully employed to probe the Anderson localization in the 
asymptotic time regime Hi]. Moreover, it is possible to 
imprint various homogeneous disorder distributions on 
the speckle potential . The restriction to short times 
would be implemented by simply switching the speckle 
potential off after the desired exposure time. Time-of- 
flight measurements then reveal the momentum distri¬ 
bution of the state. Producing an initial spatial super¬ 
position state, one may in this way observe the disorder- 
induced transition from a superposition into a mixture in 
terms of the loss of visibility of the interference pattern 
in momentum. 


VII. CONCLUSIONS. 

We developed a theory which describes the ensem¬ 
ble average dynamics of disordered quantum systems 
at short times in terms of Lindblad master equations, 
with the statistical properties of the disorder potential 
encoded in the Lindblad terms. While this effective 
evolution equation accurately captures the onset of the 
disorder-induced coherence loss of the ensemble-averaged 
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FIG. 2. (Color online) Dynamics of a Gaussian initial state 
in the ID Anderson model: Comparison of the time evolution 


pme(t) predicted by the disorder master equation (101 with 
the numerically exact dynamics pens(t), obtained by averag¬ 
ing over a finite number of realizations, (a) Density matrix of 
the evolved ensemble average state pens(t) at t = 0.2h/J with 
strong disorder, W = lOJ, and K = 200 realizations. The 
state displays a spatially homogeneous decay of the coher¬ 
ences, as predicted by the master equation, (b) This is also 
confirmed by the ratio of the two density matrices, which re¬ 
mains everywhere close to one, with local fluctnations on the 
order of a few percent, due to the finite sampling, (c) The 
short-time approximation becomes poor beyond t = 0.2fi/J, 
which can easily be seen by inspecting the ratio of the two pu¬ 
rities Pens(t) and pine(t), which starts to increasingly deviate 
from one. The purity provides a global measure of the deco¬ 
herence and is robust, in the sense that it averages out local 
fluctuations due to the finite sampling, (d) At t = 0.2h/J the 
initially pure state has lost about 45% of its purity. The latter 
continues to decrease monotonically and eventually converges 
to Pens = 0.074 beyond t = 0.5h/ J, reflecting the remaining 
coherence in the asymptotic state. 


state in the ID Anderson model, our theory is not yet 
capable to explain other disorder effects such as diffu¬ 
sive propagation or localization. However, a (translation- 
covariant) master equation which also captures the en¬ 
semble average dynamics of such disorder-induced phe¬ 
nomena must, in principle, exist. These must then 
emerge as a feature of the, in general, incoherent evo¬ 
lution of the ensemble-averaged state. Indeed, Fig. [^d) 
illustrates the monotonic decay of the averaged state’s 
purity towards that of the (localized) asymptotic state. 
The asymptotic value of the purity decreases with in¬ 
creasing disorder strength W, which reflects that the re¬ 
maining coherence in the asymptotic state is related to 
the localization length ^ oc l/W^ |22) . 

Our results represent a hrst step towards a treatment of 
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disordered quantum systems in terms of quantum master 
equations. The impact of spectral and of unitarily invari¬ 
ant disorder on the dynamics of the ensemble-averaged 
state of finite-dimensional quantum systems at arbitrary 
times t is the subject of [25] , 
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IX. APPENDIX 


In the following we evaluate the short-time disorder 
dynamics Q for two simple, yet relevant examples: a 
particle of mass m in one-dimensional, continuous space, 
subject to a random i) linear or ii) harmonic potential. 


i) In the linear-potential case, we consider a Hamil¬ 
tonian of the form Hg = fp'/2m-\-ex, i.e. the randomness 
lies in the strength of the constant force exerted on the 
particle. This describes for example experiments where 
a charged particle is exposed to a homogeneous, but not 
fully controlled electric field, i.e. the field strength varies 
from run to run. If we assume for simplicity that e = 0, 
the average Hamiltonian corresponds to the free Hamil¬ 
tonian, H = p^/2m, and for the Lindblad operators we 
obtain Lg = £x/Eq, with Eq an arbitrary energy scale 
(which is again introduced for dimensional reasons and 
irrelevant for the final result ( [I^ ). Since all Lindblad op¬ 
erators are proportional to x, we can perform the disorder 
integral in 0 and obtain the simplified master equation 


P 




hP 


t[x,[x,p]]. 


( 12 ) 


This is the well-known Caldeira-Leggett master equa¬ 
tion [HI HU, which usually emerges in an open-system 
context from a linear coupling model. The incoherent 
part of predicts an exponential decay of spatial co¬ 
herences according to (as for the derivation of the col- 
lisional decoherence localization function 0 , we ne¬ 
glect for the moment the von Neumann commutator) 

p^{x,x') = exp —^t^{x — x'Y pq{x,x'). While such 
a localization rate which grows above all bounds for 
|x — x'l —>■ c» is usually considered as unphysical in the 
open-system context, it arises here as a natural and un¬ 
avoidable consequence of the disorder average. 

ii) In the harmonic-potential example, we could in 
principle allow for both a random frequency and a ran¬ 
dom center point. We focus here on the latter and keep 
the frequency fixed, Hg = fP/2m {muPI2)(x — e)^. 
This may describe experiments where a particle is har¬ 
monically trapped, but where the trap center is sub¬ 
ject to fluctuations. In this case (again assuming that 
e = 0), the short-time dynamics 0 result, once again. 


in Caldeira-Leggett decoherence, p = —{i/h)[H,p] — 
/hp)t[x, [a;,p]], but this time with a harmonic av¬ 
erage potential H = pPI2m -I- {moPI2)SP. In this exam¬ 
ple, we can even anticipate the evolution of the ensem- 








7 


ble average beyond the short-time approximation: Since 
all random potentials share the same frequency ui, any 
initial state will at multiples of the period T = 27r/w re¬ 
cur, and in particular it will regain the purity lost in the 
early stage. On the level of the disorder master equation, 


this indicates periodic, partly negative decoherence rates 
7 e(t). While such time dependence comprising (at least 
partial) purity revivals is likely the generic pattern of the 
ensemble average dynamics. Fig. [^d) indicates that the 
Anderson model exhibits a strictly monotonic decay of 
coherences, also beyond the short-time approximation. 


